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Abstract 

We theoretically demonstrate that detectors endowed with internal gain and operated in regimes in which 
they do not necessarily behave as photon-counters, but still ensure linear input/output responses, can allow 
a self-consistent characterization of the statistics of the number of detected photons without need of know- 
ing their gain. We present experiments performed with a photo-emissive hybrid detector on a number of 
classical fields endowed with non-trivial statistics and show that the method works for both microscopic 
and mesoscopic photon numbers. The obtained detected-photon probability distributions agree with those 
expected for the photon numbers, which are also reconstructed by an independent method. 

PACS numbers: 42.50.-p (quantum optics), 42.50.Dv (nonclassical states), 42.50.Ar (photon statistics and coherence 
theory), 03.65.Wj (State reconstruction, quantum tomography) 
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INTRODUCTION 



Full characterization of quasi-monochromatic electromagnetic fields requires the knowledge of 
amplitudes and phases, which can be achieved by using homodyning detection techniques. These 
techniques, which belong to routine in radio-frequency analysis, were brought to operation at op- 
tical frequencies LI] and became popular since the discovery of phenomena generating fields with 
non-classical properties (see Ref. [2] for a review). They were used not only in the case of quasi- 
monochromatic light [3], but also in the case of definitely broadband fields, that is of pulsed fields 
|4|]. However, a simple information concerning the field amplitude such as a statistical informa- 
tion, for instance of the photon-number, is often sufficient and sometime crucial to understand 
light-matter interaction phenomena [51]. 

It might be thought that measuring photon-number statistics is feasible only when very few 
photons are present in the field to be characterized [6]. This prejudice stems from the notably poor 
photon-number discriminating capability of detectors, which hardly goes beyond five detected 
photons 112, l8|] for detectors based on the photoelectric effect. Many efforts have been made to 
achieve better photon-counting capabilities by working both on the detectors and on the front-end 
optics. Multi-pixel detectors with single-photon sensitivity in each pixel, such as arrays of single- 
photon solid-state detectors [9], intensified CCD cameras |Ji3] and silicon photomultipliers 1 11], 
were (are being) shown to allow photon counting provided the light is spread across the sensitive 
area so that each pixel receives one photon at most. Photons temporally spread by either cascaded 
arrays of beam splitters or multiple fiber- loop splitters have been alternatively used in connection 
with single-photon avalanche diodes [1^]. Some of us showed that an on/off detector, that is a 
detector giving a standard current output as the response to any n > 1 photon number, allows 
determining the whole statistical distribution of photons from the experimental probability of de- 
tecting zero photons at varying the detection efficiency. The technique, based on the application 
of a maximum likelihood algorithm, has been exploited to reconstruct both classical and quantum 



states (|l3L ll4l]. However, the method loses efficacy at increasing n, nor the result improves if the 
detector has photon-number discriminating capabilities for n > 1 [ | 13f | - 

In this work we use a hybrid photo-detector (HPD), which is a photoemissive detector endowed 
with internal gain as it includes a photocathode and an avalanche diode. The pulse-height spectra 
of the output charge of this HPD exhibit separated peaks typically for m < 3 — 4 photoelectrons. 
We perform single-shot measurements on ps-pulsed classical fields producing mean charge values 
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from below to well above such m values. From the converted voltages, v, collected for thousands 
of measurements, we determine the ingredients of the Fano factor of the voltages, i.e. mean and 
variance of the P v distribution. We show that, if the detector is operated within its range of linear 
input/output response, we can retrieve P^ from the experimental t>-data without need of cali- 
brating either the HPD or the electronics that processes its output. In the cases of relatively low 



intensities, in which the P el peak is sizeable, we apply the algorithm cited above A13I1 to recon- 
struct the statistical distribution of photons and find P ph distributions in good agreement with the 
theoretical ones. It is noticeable that the method easily allows determining the statistical distribu- 
tion of n < 100 photons, that is in the mesoscopic regime, which has been scarcely explored for 
the above reasons. 

The paper is organized as follow. In Section U we develop the theory that allows the retrieval 
of the detected-photon distribution, P^ 1 , from the experimental P v . In Section HH we check the 
theoretical results with experiments performed on a number of non-trivial classical fields with our 
hybrid photo-detector. For each type of field, in Section [Till we reconstruct the distributions of the 
photon-numbers P ph from the simple knowledge of the P el peaks of the distributions retrieved 
in Section [III In the conclusive Section we present some perspectives of our method as one that 
allows determining photon-number statistics with linear non-photon-counting detectors. 

I. THEORY 

For any photo-detection experiment performed with an overall detection efficiency 77 < 1, the 
distribution of the number of electrons generated by the primary photon interaction that occurs in 
the detector, PJ l , is linked to that of the number of photons in the light under measurement, P n ph , 
by 
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If we limit our analysis to the first two moments of the distributions, the link between the statistics 



of photons and that of the detected photons is given by lllql 
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rjn ; = 7] 2 a {2) +rj(l-r])n, (2) 



in which n = J2^=o n Pn ph * s me mean value and = Yl^=o( n ~ n) 2 P n ph is the variance. The 
same notation is adopted for fh and 



r (2) 
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In principle, P^ 1 ^ P n ph , that is the distribution of the photoelectrons is not necessarily the 
same as that of the photons. As a matter of fact, all the fields studied in this work are endowed 
with a photon-number distribution that is invariant under the convolution in Eq. ([]])• We will exploit 
this property to compare the experimental results with the corresponding theoretical distributions. 

For photoemissive detectors, in which the primary detection process is the photoelectric effect, 
m is the number of photoelectrons emitted by the photocathode, which normally undergo multipli- 
cation before reaching the anode. For the linearity of the response, both the amplification internal 
to the detector (gain) and that of the electronics that processes its output, must be independent of 
m. In the following we also adopt the strongly simplifying hypothesis that the spread of the single 
photoelectron peak in the final electronic output is negligible as compared to its mean value. Un- 
der such an hypothesis, which is too weak to ensure photoelectron counting capability, the relation 
linking the statistics of the photoelectrons to that of the final voltage outputs, v, is 

7 

being 7 a conversion coefficient describing the overall m-to-t> conversion process. The distribution 



in Eq. © is the experimental P v . Its first two moments are given by Ill6ll 

v = 1 m ■ a™ = , (4) 

where the symbols are defined as in Eq. ©. 

Apparently Eq. © allows retrieving P^ 1 from the experimental P v only if 7 is known. The 
problem is that to obtain the 7-value one should repeat the measurement of P v for the same state 
of light at different mean numbers of photons n by using a detector with photoelectron counting 
capability. We show that we can avoid such direct calibrations and determine the 7-value by other 
means. In fact, if we calculate the Fano factor for the output voltages, F v , by using Eq. © and 
Eq. © we find 

_ c4 2) 7 2 [i] 2 a^ + 77(1 - 17)71] 
r v = — = z 

v 77771 

= 777^ + 7(1 - 77) = 777Q + 7 , (5) 

(2) 

where we have rewritten the Fano-factor for the photons, F = a^[n)/n, as F = 1 + Q by 



inserting the Mandel Q-factor II 1711 . By multiplying and dividing Eq. © by n and reusing the 
above equations we get 

Q 

F v = ^-v + 1} (6) 

n 



in which Q/n = (cr^(n) — n)/n 2 depends only on the specific state under measurement. Thus 
Eq. © shows a linear dependence of the F v -factor on v, with Q/n as the proportionality coeffi- 
cient, which is zero for Poissonian light, positive for classical super-Poissonian light and negative 
for nonclassical sub-Poissonian light. This linearity can be verified by repeatedly measuring a field 
upon inserting filters with different transmittance in front of the detector. In fact the insertion of 
filters amounts to varying v, but leaves unaltered the quantity Q/n. In practice, as the expression 
in Eq. © is general, we expect that, for any given statistics, upon plotting the experimental values 
of F v as a function of v the data points should align along a straight line with intercept 7 and slope 
Q/n. 

In order to exploit Eq. © to retrieve the P^ distributions for lights with different photon- 
number statistics, we built some "artificial" field states by mixing, at a beam splitter, different 
fields generated as described in Section lU With reference to Fig. Q] we anticipate a summary of 
the schemes we adopted to produce the field states to be measured: (A) coherent light directly from 
the laser source; (B) single-mode thermal light obtained by passing the coherent field through a 
diffuser, namely a rotating ground glass plate (D in Fig. [Q) and selecting a single speckle by a 
pin-hole (PH in Fig.Q]); (Q multi-mode thermal light obtained as in (B), where multiple speckles 
where selected by a wider pin-hole and focused by lens L; (D) phase-averaged displaced coherent 
light obtained by mixing two coherent fields at beam splitter BS upon averaging the relative phase 
between them by using a piezoelectric movement (Pz in Fig.Q]); (E) displaced single-mode thermal 
light obtained by mixing the state generated in (B) with a coherent field; (F) displaced multi-mode 
thermal light (same as (E) with state (C) instead of (B)). In cases (A), (B) and (C), the neutral- 
density filter F in Fig. [His substituted by a beam stop. In cases (D), (E) and (F), F is inserted on the 
coherent beam. As states (D), (E) and (F) are characterized by non-trivial statistics, reconstructing 
P^ 1 for such states would be a rather convincing test of the goodness of our method. 

In the following we present the expressions of the statistical distributions, P n , of the photons 
for the above states, for which we calculate the slope coefficients Q/n of F v as a function of v 
in Eq. ©. All mean photon numbers appearing in the following calculations are taken at the BS 
output. 
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FIG. 1: (Color online) Schematic diagram of the experimental setup. HPD, hybrid photo-detector; SGI, 
synchronous gated integrator; F, neutral density filter; BS, 50% beam splitter; L, focusing lens; P, polarizer; 
Pz, piezoelectric movement; D, diffuser (rotating ground glass plate); PH, pin-hole. Components in dotted 
boxes are inserted/activated when necessary. 



A. Coherent 

The photon-number distribution of the coherent field is given by the Poisson distribution 
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for which n = = \a\ 2 , so that, according to Eq. ©: 



F v = 7 



(7) 



(8) 



independently of the mean value v. 



B. Single-mode thermal 



The photon-number distribution of a single-mode thermal field is given by 
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for which n = nth an d o\ = n t h {nth + 1)- In this case, Eq. © gives 

F v = v + 7 ■ 



(9) 



(10) 



C. Multi-mode thermal 



For jx independent thermal modes equally populated by N t h / [i mean photons Q6Q we have the 
distribution 

(n + LL-iy. 
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for which n = N th and a 2 = N th (N th / fi + 1). Equation © gives 

F„ = - + 7 . (12) 
H 

Note that the slope coefficient 1/// has a value lying between and 1 . 
D. Phase-averaged displaced coherent 

The state is the superimposition of a coherent state with |ai| 2 mean photons and a phase- 
averaged state with |a 2 | 2 mean photons at the same frequency. The expected photon-number 



distribution for such a state is given by [|14|1 

A n A A ( n\ (-l) k fB\ k r(l/2 + fc/2) T (1/2) rr , , r , 5 

* = E^J^U) T^r - f - [ {1/2 + */ 2 ' ■ w 2 - 1 + • B ' 

(13) 

where A = |ai| 2 + |ck 2 | 2 , £> = 2|ai||a 2 | and iF 2 (a, b, z) is the generalized hypergeometric func- 
tion. The distribution in Eq. (fT3l) has mean n = |«i| 2 + |ct2 1 2 and variance = n (Kn + 1), 
withiT = 2\a 1 \ 2 \a 2 \ 2 /(\ai\ 2 + \a 2 \ 2 ) 2 . We thus have: 

F v = Kv + 1} (14) 

in which the value of the slope, K, lies between and 1/2. This fact makes the dependence on v 
in Eq.(fT4l) not distinguishable from that in Eq. (fT2l) when ji > 2. 



E. Displaced thermal 

The state is the superimposition of a thermal state having n th mean photons and a coherent state 
with | ct | 2 mean photons. The photon-number distribution is 
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■n- +1 -exp — )L n [ —— J , (15) 

{n th + l) n+v \ n th + lj V n th (n th + l)J 

in which L n is the Laguerre polynomial for 7 = 0. In this case we find n = n th + \a\ 2 , 
<Jn^ = n th + I ct 1 2 + n th (n th + 2|a| 2 ) and finally: 

F v = K l v + 1 , (16) 

where K x = n th (n th + 2|a| 2 ) / [n th + |a| 2 ) 2 . Note that the value of K\ lies between and 1 on 
varying \a\ 2 . 



F. Displaced multi-mode thermal 



The state is the superimposition of the state in subsection (C) with a coherent state with \a\ 
mean photons. The photon-number distribution is given by 

p - jWg (_ H 2 \ r ,-t (_ H 2 \ nTk 

for which we have n = iV^ + ,u|a| 2 and = N th + /i|a| 2 + N th (N th //j, + 2|a| 2 ). Thus we find: 

F v = K 2 v + j, (18) 

where K 2 = K'Jfj,, being K[ defined as K[ = N th / /i {N th /fjL + 2|a| 2 ) / {N th / /i + \a\ 2 ) 2 in anal- 
ogy with case (E). In this case, on varying \a\ 2 , the slope K 2 takes values between and 

Note that the slopes K and K\ in the above cases (D) and (E), respectively, can be simply 
expressed in terms of the ratio of the intensities of the mixed fields, while for K 2 the knowledge 
of the number of modes fj, is also necessary. We finally remind that all the above distributions 
are identical to those of photoelectrons, upon rescaling the parameters by the overall detection 
efficiency, 77. 



II. EXPERIMENTAL RESULTS AND DISCUSSION 



All the measurements presented here were performed on picosecond-pulsed light fields at 523 
nm wavelength. The light source was the second-harmonic output of a Nd:YLF mode-locked laser 
amplified at 5 kHz (High Q Laser Production, Austria) producing pulses of ~5.4 ps duration. We 
used a multimode fiber with 600 fim core diameter to deliver the light to the HPD detector (H8236- 
40, Hamamatsu, Japan, maximum quantum efficiency of the photocathode: 0.40 at ~550 nm). The 
maximum overall detection efficiency, calculated by including the losses of the collection optics, 
was rjmax = 0.29. The detector was strictly operated within its range of linear response. Its output 
current pulses were suitably gate-integrated by SR250 modules (Stanford Research Systems, CA) 
and sampled to produce a voltage, v, which was digitized and recorded at each shot. 
As sketched in Fig.[Q we used a 50% non-polarizing beam splitter, BS, to mix the fields in cases 
(D), (E) and (F). In case (D) the phase randomization from shot to shot was obtained by changing 
the path of one of the two fields at a frequency of ~ 100 Hz with a piezoelectric device, Pz, covering 
a travel length of 1 .28 /mi. Pseudo-thermal light was obtained by means of a diffuser, D. Less than 
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one speckle matched, in cases (B) and (E), the head of the fiber delivering the light to the detector. 
The insertion of lens L together with the pin-hole PH provided the coupling to the fiber of multiple 
speckles in cases (C) and (F). 

To obtain the F„-value as a function of v, we inserted a polarizer, P, as shown in Fig. [TJ This 
amounts to changing the overall detection efficiency r\. In each of the cases (A) to (F), we recorded 
at least 30000 w-values for each angle of the polarizer and repeated the measurements upon blank- 
ing either of the BS inputs. We also measured one set of data in the absence of light, whose mean 
value was adopted to set the origin of the ?;-values (i.e. v = 0) of the data collected for each light 
field. Note that, using a single detector and a single laser and not moving any optical component 
throughout all the experiments make 77 virtually constant if not for the changes in the polarizer 
rotation angle. 

For each set of data we calculated mean value, v , and variance, ai, , thus obtaining the i^-value 
corresponding to each P v distribution, upon correcting for the variance of the data set mea- 
sured in the dark. 

The insets in Fig. |2] display the plots of F v as a function of v we obtained in cases (A) to (F). 
In agreement with Eq. © the experimental data are well fitted by straight lines, whose intercepts 
at v = give 7. We find the following values of 7, in Volt: 0.21 (A), 0.20 (B), 0.21 (C), 0.21 
(D), 0.18 (E) and 0.19 (F). The values of the slopes of the best fitting straight lines, which are 
in (A), 0.980 in (B), 0.189 in (C), 0.491 in (D), 0.863 in (E) and 0.054 in (F), can be interpreted 
by using the results in the Section HI Note that in cases (A) and (B), which correspond to coherent 
and single-mode thermal lights, respectively, we virtually find the slopes expected, that is zero and 
one (see Eq. © and Eq. (flOl) ). 

Once 7 is determined, we can reconstruct the photoelectron statistics, P„f . We convert voltages 
v into number of photoelectrons m by dividing the v -values by the value of 7 according to Eq. ©. 
We then rebin the obtained distributions in unitary bins to find . In the main panels of Fig. [2] 
the bar plots are the reconstructed photoelectron distributions for some of the w-data sets used to 
obtain the calibrations in the insets and the white-bar plots correspond to the cases of the maximum 
v-values. For the last ones we have the following mean values: fa = 1.95 in (A), 0.60 in (B), 2.07 
in (C), 2.41 in (D), 1.15 in (E) and 4.25 in (F). The bar plots in grey/black correspond to different 
(lower) m-values (data not shown) and equal ratios between the mixed fields. 

In the case of the coherent field (A), knowing m is enough to calculate the theoretical distribu- 
tions of photoelectrons by using Eq. ([?]), which are displayed as full lines in Fig. [2] (A). To assess 
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the quality of the reconstruction we calculate the value of the fidelity: 

oo 
3=0 

For all the reconstructions in the figure we find / = 0.9996 ± 0.0004. In particular the reconstruc- 
tion of the white-bar plot having the above-mentioned m yields / = 0.9996. 

The case of single-mode thermal light (B), is similar. By using the experimental m-values into 
Eq. © we obtain the curves displayed in Fig. [2](B) and / = 0.9993 ± 0.0004. In particular 
/ = 0.9992 in the case of maximum intensity. 

Cases (A) and (B) are the only ones for which the slope, Q/n, of F v as a function of v identifies 
the statistics of the light. The proximity of the slopes of the fitting straight lines to and 1, 
respectively, makes the high values of / a result to be expected. 

In case (C), in which the field is the superimposition of fi equally populated thermal modes, 
in order to calculate the theoretical distributions by means of Eq. (fTTI) we need the additional 
information on the number of modes. However, this information is provided by the slope of F v . 
For the data displayed in Fig. |2](C) we have fj, = 5.3. The values of the fidelity are in the range 
/ = 0.999 ± 0.001 and, in particular, / = 0.9993 for the most intense field. 

Calculating the theoretical distributions for the mixed field states in cases (D), (E) and (F), 
requires the further knowledge of the relative populations of the two fields. We get this information 
from the measurements that we repeated upon blanking either of the BS inputs. 

In case (D) we determined the ratio R = Vi/ifa ~ 0.927, which remains the same for all 
the data in Fig. |2] (D), from the independent measurements of each of the coherent beams. By 
using this ratio and the experimental mean values in Eq. (fl~3l) we calculated the theoretical curves 
displayed as full lines in Fig.|2](D) and found fidelity values in the range / = 0.9995 ± 0.0003 and 
/ = 0.9992 for the most intense field. By writing the slope K in Eq. CE1) as K = 2R/(R + l) 2 
and using the above value of R, we find K = 0.494, to be compared with the slope of the fitting 
straight line in the inset (K = 0.491). 

In case (E) we measured a ratio R = v th /v a = 1.535 between the thermal field and the coherent 
displacement field. The theoretical curves in Fig.[2](E) were calculated by inserting this ratio and 
the experimental mean values in Eq. (fl~5l) . The corresponding fidelity values are in the range 
/ = 0.998 ± 0.002, being / = 0.9992 for the most intense field. The slope K x in Eq. (O, which 
turns out to be K x = R(R + 2)/(R + l) 2 , becomes K\ = 0.844 by using the above value of 
R. Note that the slope of the fitting straight line is K\ = 0.863. These results are satisfactory 
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even if the experimental data points representing F v are rather noisy (see inset in Fig. [21(E)). Note 
that also the corresponding data in the inset of Fig. [2] (B) are rather noisy, a fact that testifies the 
difficulty in producing a reliable single-mode thermal field. 

Finally, in case (F), to calculate the theoretical distributions we need to know both the ratio 
R = Vth/v a between the multi-mode thermal field and the coherent displacement field and the 
number of modes \i. The theoretical curves in Fig.[2](F) were calculated by using, in Eq. (fTTT) . the 
experimental value R = 2.446 and setting fj, — 8. This value of ji optimizes the fidelity in the case 
of the most intense field. As a check of the choice of the parameter /i, we compare the value of the 
slope K 2 = R/fi(R + 2fi)/(R + fi) 2 in Eq. £[8]) as obtained by using R and fi, K 2 = 0.052, with 
that evaluated from the fit of the Fano factor, K 2 = 0.054: the result is satisfactory, as confirmed 
by the fidelity values, which are in the range / = 0.9992 ± 0.0003, being / = 0.9990 that of the 
most intense field. 



1.0- 
0.8- 
0.6- 
0. 0.4 
0.2- 
0.0- 









_ 0.3 


















0.0 0.2 0.4 

»(Volt) 














5 






0123456789 10 11 12 1 2 3 4 5 6 7 8 

m m 




0123456789 1 2 3 4 5 6 7 0123456789 

1.0 -t 



01 23456789 10 11 12 



FIG. 2: Insets: Fano factor, F v , as a function of v, for the different light states. Bars: reconstructed pho- 
toelectron distributions, P^ 1 , for some of the data sets used to calculate the Fano factor. Lines: theoretical 
curves. 
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III. RECONSTRUCTION OF THE PHOTON STATISTICS 



It has been shown that suitably elaborating the probability of zero detected photons, P el , at 
varying the overall detection efficiency, 77, can lead to the reconstruction of the photon-number 
distribution P n ph fl 1 3M . The data needed to pursuit this approach can be collected by using a detector 
whose output simply allows discriminating the case in which no photon is detected (m = 0) from 
any case in which m > 1. Obviously our HPD has such a feature. 

We note that the data points in each of the insets in Fig. [2] were obtained from measurements 
performed at varying i] up to r\ max . We thus revisited our experimental data to the maximum likeli- 
hood method to recover the photon-number statistics as described in II 1 3M - To this aim, we took the 
values of P el from all the experimental P^ distributions and assigned the corresponding values 
of 77 as 77 = r] max v/v max . The data were then processed by the maximum likelihood algorithm 
giving an iterative solution for P ph corresponding to r] max 111311 - I n the insets of Fig. [3] we plot P el 
as a function of 77 and, as full line, the theoretical behavior expected for the field considered. We 
point out that the almost ideal dependence of our P e '-values on 77, which can be observed in the 
insets, is not enough to guarantee high fidelity reconstruction II 1 3M . 

In the main panels of Fig. [3] we show, as bars, the reconstructed P n ph distributions that cor- 
respond to the photoelectron distributions displayed as white bars in Fig. [2l The reconstruction 
algorithm was stopped after a number of iterations, which was different from state to state, giving 
the reconstructed P n ph with n- value best matching the expected one, n max = v max / ( , yr] rnax ). The 
iterations were: 1000 in (A), 1500 in (B), 1000 in (C), 48000 (D), 1500 in (E) and 10000 in (F). 

In each panel in Fig. [3] the dots show the plot of the corresponding theoretical P n distribution 
(see Section IJ) as calculated by using the experimental values of n max , R and /1 (in cases (C) and 
(F)). The /-values evaluated in analogy with Eq. (fT9l for the reconstructions in Fig.[3]are: 0.9988 
(A), 0.9993 (B), 0.9960 (C), 0.9975 (D), 0.9973 (E) and 0.9950 (F). It may be noted that these 
values are still rather close to one, but lower than those for the corresponding reconstructed P^ 1 
distributions (see white-bar plots in Fig. [2]). 

In particular, we note that the reconstruction in cases (C), (D) and (F) is not as good as in the 
others. This is mainly due to the fact that cases (C), (D) and (F) represent the cases with the highest 
mean number of photons for which we can have fewer significant (i.e. non-zero) values of P el to 
be used in the reconstruction algorithm. In addition, in case (D) the low quality of the reconstruc- 
tion confirms the recognized difficulty of the method in recovering multi-peaked distributions [15]. 
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An improvement in this direction could be given by an increase in the overall quantum efficiency. 
Moreover, in cases (C) and (F), additional problems arise from the preparation of the multi-mode 
thermal state to be measured, which, as described above, was generated by means of a rotating 
ground glass plate, selected by an aperture and focused by a lens. Any instability affecting the col- 
lection of the light by the fiber is more detrimental in this case. On the contrary, the multi-mode 
thermal state produced by a spontaneous down-conversion process we presented in ft 1 311 , which 
was more robust, could be better reconstructed. We finally note that the self-consistent method 
for reconstructing the photoelectron distributions we presented in this paper is not as much 
sensitive to noise and works very well also in cases (C) and (F). 




FIG. 3: Bars: reconstructed photon-number distributions, P n , for the different light states. The displayed 
reconstructions correspond to the measurements performed with rj = rj max , whose P^j are plotted as white 
bars in Fig. |2] Dots: theoretical curves, P n (the connecting lines are a guide for the eye). Insets: values of 
P el as a function of rj (dots) and theoretical behavior expected for each field (lines). 
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CONCLUSIONS 



We have demonstrated that it is possible to implement a procedure to recover the distribution 
of detected photons that avoids calibration both of the photo-detector and of the electronics pro- 
cessing its output. The data to be treated are the voltages, into which the amplified charge of 
the photoelectrons are converted, which are the most straightforward results of this kind of mea- 
surements. The procedure employs the evaluation of the Fano factor of the voltages measured at 
different values of the overall detection efficiency. It has been satisfactorily tested on classical 
states endowed with non-trivial statistics. In all the cases we achieved high-fidelity reconstruc- 
tions of the photoelectron statistics and found results in good agreement with the photon-number 
statistics, independently reconstructed. Works are in progress to apply the method presented here 
to nonclassical fields, such as squeezed states and conditionally prepared states. We expect to be 
able to characterize the nonclassicality of the state by intensity measurements. 
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